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We consider a massive relativistic Bose gas with N complex scalars at finite density. At zero 
temperature, we calculate the pressure, charge density and the speed of sound in the one-loop 
approximation. In the nonrelativistic limit, we obtain the classic results for the dilute Bose gas. 
We also discuss finite-temperature effects. In particular, we consider the problem of calculating 
the critical temperature for Bose-Einstein condensation. Dimensional reduction and effective-field- 
theory methods are used to perturbatively calculate the effects of the nonstatic Matsubara modes. 
Calculations of T c in the effective 3d theory require nonperturbative methods. Using the Monte Carlo 
simulations of X. Sun [Phys. Rev. E67, 066702 (2003)] and the seven-loop variational perturbation 
theory (VPT) calculations of B. Kastening [Phys. Rev. A 70, 043621 (2004)], we obtain T c for 
N = 2 to second order in the interaction. 

PACS numbers: 



I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) of 
trapped alkali atoms almost ten years ago has created an enor- 
mous interest in the properties of the weakly interacting Bose 
gas QLQiQ- The temperature at which these systems Bose 
condense is of the order 10 -5 Kelvin, which is many orders 
of magnitude higher than the condensation temperature and 
the onset of superfluidity in 4 He (2.17K). Both trapped alkali 
gases and 4 He are examples of nonrelativistic systems. 

BEC in relativistic Bose systems typically takes place in 
matter under extreme conditions. For example, kaons may 
condense in the color-flavor locked phase of high-density 
QCD H H 0. This phase is a superconducting phase of 
QCD which arises from an instability of the Fermi surface; In 
analogy with ordinary BCS-theory, a weak attraction among 
the quarks in one or more channels results in the formation 
of Cooper pairs and the spontaneous breakdown of the color 
symmetry of QCD 0, 0, HI- Such a phase may be found 
in the interior of compact stars if the density is sufficiently 
high. The linear SU(2)l X 5'L r (2)_R-symmetric sigma model 
at finite chemical potential \i for the hypercharge is used as 
a toy model for the description of kaon condensation in the 
color-flavor locked phase of QCD [Tol ITU. 

Calculation of the critical temperature or the critical den- 
sity for Bose-Einstein condensation of an ideal Bose gas 
has been a standard text-book calculation for a number of 
years |T^|. The charge density n of excited bosons as a func- 
tion of temperature is given by 
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where w = \/p 2 + m 2 , (3 — 1/T and il is the chemical poten- 
tial. We have set h = ks = 1. Bose-Einstein condensation 
takes place when the chemical potential is equal to the mass 
m of the bosons. At that temperature, all the charge can no 
longer be accomodated in the excited states and it condenses 
into the ground state. Generally, T c can only be calculated 
numerically, but in various limits analytic results can be ob- 
tained. For example, in the nonrelativistic (NR) limit, the 



critical density n c calculated from Eq. Q is 
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Inverting this equation, one obtains the well-known result 
T c = 2ix/m [n/C (§)] . Similarly, for m -> 0, one finds 
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Thus, in the ultrarelativistic limit m — 0, T c is infinite, or 
equivalently, the critical charge density is zero. All charge 
resides in the ground state irrespective of the temperature. 

In the context of nonrelativistic field theory, the problem 
of calculating the transition temperature for Bose-Einstein 
condensation with a weak interaction has a very long history 
and conflicting results have appeared in the literature Q. A 
first-order perturbative calculation gives no correction to the 
ideal-gas result, while higher-order calculations are plagued 
with infrared divergences. This is a typical example of in- 
frared divergences that arise in the vicinity of a second-order 
phase transition. The long-distance physics in the critical re- 
gion is nonperturbative and one has to sum up an infinite set 
of diagrams to obtain a finite result. The problem was solved 
only recently by Baym et al. [l3T ] who realized that it can 
be reduced to nonperturbative calculations using a classical 
three-dimensional field theory. Once it was understood how 
to organize the problem, the calculation of T c has later been 
carried out using several methods. These include 1/N tech- 
niques Il4l. lip . lattice simulations [TrUTH ]. the linear delta ex- 
pansion |l8lll9| , variational perturbation theory [20L I2H |22^| , 
and renormalization group methods |23t 12 If . In particular, 
the use of effective field theory methods to obtain an effec- 
tive three-dimensional field theory combined with high preci- 
sion lattice calculations has settled the issue in a very elegant 
way [IMl25|. 

The 0(2)-symmetric relativistic Bose gas at finite temper- 
ature and chemical potential has been studied in detail by 
Benson, Bernstein, and Dodelson [2(1 l^ . while the prob- 
lem of calculating T c was addressed in a paper by Bedingham 
and Evans employing the linear delta expansion |2S||. In the 



present paper, we examine the relativistic Bose gas in a more 
general setting where we consider N coupled scalars. For 
N = 1, it reduces to the standard case, while for N — 2, 
the system is isomorphic to the linear SU(2)l X SU(2)r- 
symmetric sigma model and hence relevant for kaon conden- 
sation in stars. 



stituting Eq. (J^J into Eq. 0, the action can be written as 
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II. PERTURBATION THEORY 

In this section, we briefly discuss the perturbative frame- 
work for a massive Bose gas with N charged scalars at finite 
chemical potential fi. The action is 
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where "1> = ($i, $2--., $jv) and $i is a complex scalar field. 
We first parametrize the quantum field 3>i in terms of a time- 
independent vacuum expectation value 0o and two real quan- 
tum fluctuating fields: 
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Similarly, the remaining complex fields $2, ^iv are 
parametrized in terms of 2N — 2 real fields (j>3, 02at- Sub- 
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The propagators that correspond to the free part Sf ree of the 
action are given by 
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where the dispersion relations are 
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Here the tree-level masses mi and m 2 are 
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In the minimum of the classical action, m\ = 2(fj? — m 2 ) and 
m 2 = 0, and so the dispersion relations reduce to 



ui±(p) = \Jp 2 + Sfi 2 — m 2 ± ^/ (3/7 2 — m 2 ) 2 + 4fi 2 p 2 , 

(16) 

W2± (p) = vV +^ 2 ± A* ■ 



From these equations, we see that there are two massless 
modes that in the long-wavelength behave as 
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The other excitations wf(p) and (p) are gapped with gaps 
Ai = y/2(3fi 2 - m 2 ) and A 2 = 2/z. In the case N = 2, the 
gapless particles wjf and w^T carry the quantum numbers of 
(17) K+ and while the massive modes a;* and u> 2 carry those 
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of if- and K° [I3IT1. Note that there are only iV massless 
modes despite the fact that the potential has 2N — 1 flat di- 
rections which also is the number of broken generators. This 
is in agreement with the counting rule derived by Nielsen and 
Chadba [2fl|. which states that the modes with a quadratic 
dispersion relation must be counted twice. Secondly, due to 
the quadratic dispersion relation for small p, the Landau cri- 
terion [3(J for superfluidity can never be satisfied, except for 
N = 1. Thus despite the presence of a Bose condensate, the 
system is not a superfluid. 



III. ZERO TEMPERATURE 

In this section, we apply perturbation theory at zero tem- 
perature to calculate the pressure, charge density, and the 
speed of sound in the one-loop approximation. 

The partition function Z is given by the path integral 



Z = 



where the action S is given by |0J. The pressure V is 
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where VT is the space-time volume of the system. The charge 
density can be found by differenting the pressure with respect 
to u: 
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The speed of sound c is given by the coefficient of ui-(p) 
as p — > 0. Corrections to the tree- level result can be found 
by calculating the dispersion relation in the long-wavelength 
limit including the self-energy function Hi-(ui,p). It can also 
be derived once we know the charge density and is given by 
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The chemical potential measures the amount of energy needed 
to add a particle to the system, and in the nonrelativistic 
limit, we introduce the nonrelativstic chemical potential /xnr 
by a = m + /iNR- In the NR limit, Eq. 123H is therefore 
replaced by 
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The mean-field pressure Vo is found by evaluating minus 
the classical thermodynamic potential £lo(fi,(j>o) at the mini- 
mum of the classical action So: 



W = 7X(^-™ 2 ) 2 • (25) 
The one-loop contribution to the effective potential is 

Qi(p,<h) = \i J J [lndetDi(w,p) 

+(JV-l)]ndetD a (a;,p)] +Aim 2 ^ 
+AiA^o + Ai£ , (26) 



where Aim 2 , AiA, and Aif are the one-loop mass countert- 
erm, coupling constant counterterm, and vacuum countert- 
erm, respectively. After integrating over the energy u), we 
obtain 
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The integral involving u)2± can be calculated analytically in 
dimensional regularization, but the integral of a>i± cannot. 
In order to extract the divergences analytically, we make sub- 
tractions in the integrand that render the integral finite in 
d = 3 dimensions and then extract the poles in d — 3 from 
the subtracted integrals. The substraction term Q EU b should 
not introduce any infrared divergences. Our choice for the 
subtracted integral is 
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The first two terms in f2 su b vanish identically in dimensional 
regularization since there is no mass scale in the integrand. 
The last term is given in Eq. JA.3I . The one- loop thermody- 
namics potential can then be written as 
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where L = In ) and g is a function of the ratio m/u that 
must be evaluated numerically: 



g{m/u) = - I wi± 



(30) 



The counterterms necessary to cancel the poles in e are |2 
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The one-loop contribution to the pressure Vi is given by — Q.\ 
evaluated at the classical minimum. After renormalization, 
the pressure through one loop reduces to 
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We next consider the NR limit of the pressure. In this limit, 
Mnr -C m. Moreover, the kinetic energy is much smaller than 
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m and so we can expand physical quantities in powers of the 
dimensionless quantities /^nr/wi and k 2 /2m 2 . This yields 
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The other quasiparticle excitations have uj\+ = 0J2+ = 2m and 
so their contribution can be neglected. In NR field theory, it is 
customary to set 2m = 1 and we will do so in the remainder of 
this section. Introducing the scattering length a = A/8nm = 
A/4-7T, the pressure becomes 
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Note that u>2- does not contribute to the pressure since there 
is there is no scale in the integral and so it set to zero in 
dimensional regularization. Using a simple ultraviolet cutoff 
A to regulate the integral, the divergence would be cancelled 
by a vacuum counterterm . Using Eq. IA.4I . we obtain 
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B. Charge density and the speed of sound 

We next consider the speed of sound, which is given by 
Eq. I|23|l . For simplicity we consider only the ultrarelativistic 
and nonrelativistic limits. In these limits, the tree-level results 
are c = l/y/3 and c = respectively. 

The charge density can be calculated using Eqs. g2J 
and <5H : 

n = A~ { 1 + (4^ [ 4(7V + 4)i + 2N ~ 7 + 647r2 »(°)] 

(39) 

Inverting 1391 and using 123H . we obtain the speed of sound 
due to interactions in the medium: 
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The sign of the correction is determined by the beta-function. 

In the NR limit, the charge density follows from Eqs. 12211 
and GB}: 
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The speed of sound then becomes 
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where we have eliminated fJ,\ rm NR in favor of the density n. 
Note that the expansio n pa rameter in the NR limit is the di- 
mensionless quantity \fncfi which is referred to as the the gas 
parameter. This result 1)420 was first derived by Beliaev [32| . 
who calculated the leading corrections to the dispersion rela- 
tion llMlil in the low-momentum limit 2 



IV. FINITE TEMPERATURE 

We next discuss the behavior of the system defined by 
Eq. I^J & t finite temperature. 



A. Low-temperature effects 

We first consider the thermal corrections to the pressure at 
temperatures T much lower than the chemical potential p. In 
this regime, the thermodynamics is dominated by the massless 
modes. We can then approximate the dispersion relations by 
uji (p) and u>2 (p) by their low-momentum limits 1181 and 1191 . 
The pressure in the one-loop approximation is 
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Omitting the contribution from the massive modes, neglect- 
ing the zero-temperature pieces, and using lA.14ll - lA.16ll . we 
obtain 
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Similarly, in the nonrelativistic limit, one finds 
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where we again have set 2m = 1. For N = 1, this reduces to 
the old result of Lee and Yang |3^] . 



B. Dimensional reduction 

Effective field theory methods can conveniently be used to 
organize the calcuation of physical quantities whenever two 
or more momentum scales are well separated. The conden- 
sation temperature for BEC is an ideal problem for applying 
effective field theory. At finite temperature there are two 
characteristic scales in the system. The first is the correla- 
tion length which is associated with the effective chemeical 
potential. Since the phase transition is second order, the cor- 
relation length becomes infinite at T c . The second scale is 



1 In the NR limit, all the divergences at the one-loop level are 
power divergences and hence the parameters require no rcnor- 
malization if one uses dimensional regularization. See Ref. Q 
for a thorough discussion. 



2 In the orginal derivation, Beliaev expressed his result in terms of 
the condensate density no, which is different from the total den- 
sity due to the depletion of the condensate caused by quantum 
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associated with the nonzero Matsubara modes and is of order 
T. For distances much larger than the inverse temperature 
and for temperatures sufficiently close to the critical temper- 
ature, so that the effective chemical potential is much smaller 
than the temperature, the nonstatic Matsubara modes decou- 
ple and the long-distance physics can be described in terms 
of an effective three-dimensional field theory for the zeroth 
Matsubara mode. The action for this effective theory is 
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where the dots indicate higher-order operators. / is a referred 
to as the coefficient of the unit operator and represents the 
contribution to the free energy density from the nonzero Mat- 
subara modes. The parameters in Eq. il Kit are functions of 
T and the coefficients of the underlying theory Q, and are 
renormalized due to their coupling to the nonstatic Matsubara 
frequencies. These parameters can be determined by integrat- 
ing out nonstatic modes explicitly as done by Bedingham and 
Evans in Ref. |28|. but perhaps a more streamlined way of 
calculating them is by matching static Green's functions |31| . 
For the purpose of matching, the chemical potential can for- 
mally be treated as a perturbation on the same footing as the 
quartic coupling. In the effective theory m 2 is also treated as 
a perturbation. The matching is carried out by calculating 
Green's functions perturbatively in the two theories and de- 
mand they be the same for distances R much larger than 1 JT . 
This way of calculating static correlators introduces infrared 
divergences at an intermediate stage and must be regular- 
ized. Note thate these infrared divergences that appear are 
the same in the two theories and hence they cancel in the 
matching procedure. We use dimensional regularization as 
discussed in the appendix. 
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FIG. 1: One -and two-loop vacuum diagrams for /. A 
dot indicates an insertion of the operator — (/i 2 — m 2 ) + 

In Fig. Q we show the vacuum diagrams through two loops 
in the full theory and the expression is 
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where we have omitted terms of order m 4 , fi 6 etc. The sign 
w is reminder that we are neglecting infrared physics which 



will be taken care of by the effective theory. Loop correction 
to the free energy in the 3d theory vanish since there is no 
momentum scale in the loop integrals . Thus / is directly 
given by l|47|l . 
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FIG. 2: One -and two-loop Feynman diagrams for the self- 
energy II(po,p)- A dot indicates an insertion of the operator 



The mass parameter is found by matching the two-point 
functions in the two theories at zero external momenta. The 
diagrams that contribute to the self-energy function in the 
full theory through two loops are shown in Fig. |21 The self- 
energy in the effective theory vanishes for the same reason 
as did the loop corrections to the free energy. This implies 
that the mass parameter m| is given directly by evaluating 
the Feynman diagrams in Fig. Q We obtain 
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where we again have neglected terms of higher-order terms 
Z\ is the renormalization constant for the coupling A: 
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After renormalization, the mass term reduces to 
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3 Recall that m 2 for the purpose of matching is treated as a per- 
turbation and thus the propagators are massless. 
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Note that the mass parameter has a UV divergence after 
renormalization. The remaining divergence is exactly the one 
arising at the two-loop level in the 3d effective theory (see 
also Sec. llVCt . 

Finally, we need A3 at the tree level. By comparing the 
coefficients of the operator ($^<1?J in the two theories and 
taking the different normalization of the fields into account, 
one finds 

A3 = AT . (52) 



C. Critical temperature 

After having determined the parameters in the effective 
theory, the strategy for calculating T c is as follows. First we 
determine the critical chemical potential [i c as a function of 
temperature and the critical value of the m§. Then we cal- 
culate the density as a function of T and and obtain the 
critical density as a function of T by the substituing the ex- 
pression for \i c . Finally, the critical temperature is determined 
by inverting the critical density as a function of T. 

The relation between the bare mass and the renormalized 
mass is 



2 

m 3 



m 3 . 



(N- 



1)A§ 



(4tt) 



(53) 



This relation is exact due to the fact that the effective 3d 
theory I46H is superrenormalizable. Since the bare mass is 
independent of the renormalization scale, the renormalized 
mass satisfies an evolution equation. This equation relates 
the value of m| ren evaluated at two different normalization 
points A and Ao: 

mten(Ao) = mt c „(A) + {N In ^ . (54) 

Using Eq. I)53|l . Eq. do 1 1 for the chemical potential becomes 
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The charge density is given by 
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At the critical point, the renormalized mass is by dimen- 
sional analysis proportional to A3. The case N — 2 is relevant 
to kaon condensates in stars and we therefore consider this 
case in the following. In the remainder of this section, we 
also restrict ourselves to the ultrarelativistic limit. Its value 
was determined by Sun |34| using lattice simulations: 
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where the renormalization scale was chosen to be A = A3/3. 
The critical chemical potential then reduces to 
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The expectation value cannot be calculated in per- 

turbation theory due to infrared divergences. They depend 
on nonperturbative physics and can e.g. be determined us- 
ing lattice simulations or the 1/A-expansion. At the critical 
point and for N = 2, it was computed by Sun [34| using Monte 
Carlo calculations: 
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Inserting Eqs. I|59^ and iKiOt into Eq. I57|l . the critical density 
becomes 
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Inverting this equation, we obtain the critical temperature as 
a function of the density 
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The leading-order result is the usual perturbatively calculable 
high-temperature result, while the second-order term involves 
nonperturbative physics. In contrast, the first order correc- 
tion to T c in the nonrelativistic Bose gas cannot be determined 
in perturbation theory |l3j|. Note also that in accordance with 
Eq. T c becomes infinite in the absence of interactions. 
Finally, the impressive seven-loop VPT calculations of Kas- 
tening IHH^ give 0.002586(17) and -0.002796(192) for the 
quantities in Eqs. 158H and 1601 . Thus the critical tempera- 
ture is within errors in complete agreement with the lattice 
prediction. 



The quantity ( ($ t <I>) ^ is by dimensional analysis propor- 
tional to A3. Its contribution to the density is therefore third 
order in the interaction and can be omitted in a second-order 
calculation. Using Eqs. 148L I51L and 1561 . the charge den- 
sity becomes 
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V. SUMMARY 

In this paper, we have discussed the thermodynamics of 
relativistic Bose gases at zero and finite temperature. At zero 
temperature, thermodynamic quantities can be expanded in 
a loop expansion, and we calculated the pressure, charge den- 
sity and speed of sound in the one- loop approximation. In the 
nonrelativistic limit, one easily obtains the standard results 
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for the dilute Bose gas. In the critical region, perturbation 
theory breaks down due to infrared divergences. However, 
one can take advantage of the fact that there is a separation 
of scales to simplify the problem of calculating static quan- 
tities such as T c . The effects of the nonstatic modes can be 
caclulated perturbatively employing dimensional reduction, 
while the effective three-dimensional theory must be treated 
nonperturbatively. It is interesting to note that the leading 
correction to T c in the ultrarelativistic limit, is calculable in 
perturbation theory while in the NR limit it is not. The rea- 
son is simply that the chemical potential couples quadratically 
to <&t$ in the first case and linearly in the latter. 

By gauging the linear SU(2)l X 5t/(2)_R-symmetric sigma 
model in various ways, one obtains a number of interesting 
gauge theories flCl fllf . The Higgs mechanism and the Gold- 
stone mechanism are both realized in the conventional man- 
ner. The interest in these models is partly due to the fact 
that the rotational symmetry is broken as well and leads to 
a directional dependence of the dispersion relation which is 
linear for small wave vectors and roton-like for larger wave 
vectors. These models may therefore be of interest for con- 
densed matter systems such as superfluid helium. 
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APPENDIX A: FORMULAS 

Dimensional regularization can be used to regularize both 
the ultraviolet divergences and infrared divergences in three- 
dimensional integrals over momenta. The spatial dimension is 
generalized to d = 3 — 2e dimensions. Integrals are evaluated 
at a value of d for which they converge and then analytically 
continued to d = 3. We use the integration measure 
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(2tt) 3 



(A.l) 



where A is an arbitrary momentum scale. The factor (e 7 /47r) e 
is introduced so that, after minimal subtraction of the poles 
in e due to ultraviolet divergences, A coincides with the renor- 
malization scale of the MS renormalization scheme. 
We need the regularized integrals 
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In the imaginary-time formalism for thermal field theory, 
the 4- momentum P — (uj n , p) is Euclidean with P 2 — cj^+p 2 . 
The Euclidean energy po has discrete values: ui n = 2nnT for 
bosons, where n is an integer. Loop diagrams involve sums 
over uj n and integrals over p. We define the dimensionally 
regularized sum-integral by 
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The specific sum-integrals needed are 
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We also need to expand some sum-integrals about zero tem- 
perature. The phonon part of the spectrum then dominates 
the temperature-dependent part of the sum-integral. We can 
therefore approximate the dispersion relations u>i_(p) and 
u>2- (p) by their low-momentum limits I18t - 119j . and this gives 
the leading temperature correction. These are 
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